Resonant Photonic Quasicrystalline and Aperiodic Structures 



o 
o 

(N 



A.N. Poddubny.i'B L. Pilozzi,^ M.M. Voronov,^ and E.L. Ivchenko^ 

^A.F. loffe Physico- Technical Institute, 194021 St. Petersburg, Russia 
^Istituto dei Sistemi Complessi, CNR, C. P. 10, Monterotondo Stazione, Rome 1-00016, Italy 

We have theoretically studied propagation of exciton-polaritons in deterministic aperiodic 
multiple-quantum- well structures, particularly, in the Fibonacci and Thue-Morse chains. The atten- 
tion is concentrated on the structures tuned to the resonant Bragg condition with two-dimensional 
quantum-well exciton. The superradiant or photonic-quasicrystal regimes are realized in these struc- 
tures depending on the number of the wells. The developed theory based on the two-wave approxi- 
mation allows one to describe analytically the exact transfer-matrix computations for transmittance 
and reflectance spectra in the whole frequency range except for a narrow region near the exciton 
resonance. In this region the optical spectra and the exciton-polariton dispersion demonstrate scal- 
ing invariance and self-similarity which can be interpreted in terms of the "band-edge" cycle of the 
trace map, in the case of Fibonacci structures, and in terms of zero reflection frequencies, in the 
case of Thue-Morse structures. 

PACS numbers: 42.70.Qs, 61.44.Br, 71.35.-y 
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I. INTRODUCTION 

Quasicrystalline and other deterministic aperiodic 
structures are one of the modern fields in photonics 
researchjii^i^i^ Due to a long-range order such struc- 
tures can form wide band gaps in energy spectra as 
in periodic photonic crystals^ and simultaneously pos- 
sess localized states as in disordered media.- The sim- 
plest and most well-studied systems consisting of only 
two structural building blocks are Fibonacci and Thue- 
Morse chains, the former one being a quasicrystal. Pho- 
tonic crystals are known to allow a strong enhancement 
of the light-matter interaction, particularly, if the ma- 
terial system has elementary excitations and the light 
frequency is tuned to the resonant frequency of these 
excitations (the so-called resonant photonic crystals). 
In such systems the normal light waves are polaritons. 
The resonant properties of elementary excitations in 
quasiperiodic multilayered structures have been stud- 
ied for plasmons and spin waves, see the review Q, as 
well as for embedded organic dye molecules.— Whereas 
resonant periodic structures based on quantum wells 
(QWs) are widely investigated, both theoretically and 
experimentally^iiiaZiisai^ae^as^^ aperiodic long- 
range ordered multiple QWs (MQWs) have attracted at- 
tention quite recently.i2ii2SiiS^ 

In Ref . [m we have formulated the resonant Bragg con- 
dition for the Fibonacci MQWs. We have shown that 
the MQW structure tuned to this condition exhibits a 
superradiant behaviour, for a small number N of wells, 
and photonic-crystal-like behaviour, for large values of 
N. Moreover, in order to describe the light propagation 
in the infinite Fibonacci MQWs we have applied a two- 
wave approximation and derived equations for the edges 
of the two wide exciton-polariton band gaps (or pseudo- 
gaps) where the light waves are strongly evanescent. The 
fabrication and characterization of light-emitting one- 
dimensional photonic quasicrystals based on excitonic 
resonances have been reported in Refs. [2^l23l . The mea- 



sured linear and nonlinear reflectivity spectra as a func- 
tion of detuning between the incident light and Bragg 
wavelength are in good agreement with the theoretical 
calculations based on the transfer-matrix approach, in- 
cluding the existence of a structured dip in the pro- 
nounced superradiant spectral maximum. 

In this work we further develop the theory of aperi- 
odic MQWs with particular references to the Fibonacci 
and Thue-Morse sequences. In Sees. HIl and Hill we define 
the systems under study, present the results of the ex- 
act transfer-matrix computation in the superradiant and 
photonic-crystal regimes and make their general analysis. 
In Sec. IIVI we apply the two- wave approximation to de- 
rive analytical formulas for the light reflection and trans- 
mission coefficients. Comparison with the exact compu- 
tational results shows that the approximate description is 
valid in a surprisingly wide range of the light frequency uj, 
the number N of QWs, and the value of nonradiative de- 
cay rate of a two-dimensional exciton. In the close vicin- 
ity to the exciton resonance frequency loq the two-wave 
approximation is completely invalid. We show in Sec. |V] 
that in this region both the studied aperiodic structures 
demonstrate scaling invariance and self-similarity of op- 
tical properties. The main results of the paper are briefly 
summarized in Sec IVIl In Appendix the consistency of 
the two- wave approximation is questioned in terms of the 
perturbation theory going beyond this approximation. 



II. BASIC DEFINITIONS 

Here we present the definitions of the aperiodic MQW 
chains considered in this work. The structure consists of 
iV semiconductor QWs embedded in the dielectric matrix 
with the refractive index ni,. Each QW is characterized 
by the exciton resonance frequency loq, exciton radiative 
decay rate Fq and nonradiative damping F. We neglect 
the dielectric contrast assuming the background refrac- 
tive index of a QW to coincide with ni,. The center of the 
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m-th QW (to = 1 . . . N) is located at the point z = Zm, 
and the points z^ form an aperiodic lattice. Three ways 
to define a one-dimensional deterministic aperiodic lat- 
tice are based on the substitution rules^Si analytical ex- 
pression for the spacings between the lattice siteS)^ and 
the cut-and-project methodi^i^l 

We focus on the binary sequences where the interwell 
spacing takes on two values, a or b. Such structures can 
be associated with a word consisting of the letters A and 
B, where each letter stands for the corresponding barrier. 
The QW arrangement is determined by the substitutions 
acting on the segments A and B: 



A- 
B 



a{A) = AxAi 
> a{B) = B1B2 . 



• Aa+I3 
■ B^+S . 



(1) 



Each of the letters Ak and Bk in the right-hand side of 
Eq. ^ stands for A or B, a and /? denote the number of 
letters A and B in a{A), 7 and i5 are the numbers of A 
and B in <j{B), respectively The scattering properties 
of the QW sequence are described by the structure factor 



/(g)= lim f{q,N), 

N^oo 



(2) 
(3) 



Under certain conditions^liSS for the numbers a, /?, 7, 
S the structure defined by Eq. ([1]) is a quasicrystal, so 
that, in the limit N ^ 00, the structure factor ([3]) con- 
sists of (5-peaks responsible for the Bragg diffraction and 
characterized by two integers h and h' , 



h,h' — — 00 



2tt 



Ghh' = ^ I /i + — 



h' 



The parameter t in Eq. ([5]) is related by 



t = 1 



(4) 
(5) 

(6) 



with the numbers N^tN^ of the blocks B and A in 
the infinitely extending lattice. The value oi t ui ^ 
can be also expressed as t = 1 + (Ai — a) /7, where 
Ai = (/ + \/P1^4n)/2, I ^ a + 5 and n ^ (3^ - aS; 
for the quasicrystals n must be equal to ±1j^ The 
length d = {a — b)/t + b in Eq. ^ is the mean pe- 
riod of the aperiodic lattice. In the periodic case where 
a = b = d, the diffraction vectors reduce to a single- 
index set Gh = 2TTh/d with the structure- factor coeffi- 
cients \fh\ — 1. For a ^ b and irrational values of t, 
the diffraction vectors ^ fill the wavevector axis in a 
dense quasicontinuous way and the values of \fhh'\ lie in- 
side the interval (0,1). Note that, within the uncertainty 
~ {Nd)~^, the symbol (^29,0^^/ in Eq. (g]) is the Kro- 
necker delta: 52q.G^^i = 1 when 2q = Ghh' and zero when 



2q is detuned from the Bragg condition. The structure 
factor— defined without the prefactor in Eq. ^ is 
obtained by the replacement of the Kronecker symbol in 
Eq. (g]) by the functional {2n/d)5{2q - Ghh')- 

The most famous one-dimensional quasicrystal is the 
Fibonacci sequence ABAABABA . . . determined by the 
substitutions^! 



A^ AB,B^ A. 



(7) 



For the canonical Fibonacci lattice the ratios N^/Nb and 
a/b are both equal to the golden mean, r = (•\/5 + l)/2. 
The noncanonical Fibonacci structures with a/b ^ t are 
considered in Ref. [g^ and are beyond the scope of this 
paper. 

The substitution rule ([7]) can be generalized in many 
ways to provide other typ es of ID quasicrystals. It has 
been proved in Refs. |24| . [32 that any binary ID qua- 
sicrystal can be obtained by substitutions composed of 
different elementary inflations, e.g., 



A~^A''B, B^A, 



1,2. 



(8) 



For arbitrary values of a, /3, 7 and S the structure 
defined by Eq. ([T]) does not form a quasicrystal. For 
example, the substitution 



A^ AB,B~^ BA 



(9) 



defines the Thue-Morse lattice ABBABAAB . . . with a 
singular continuous structure factor and the mean period 
d = {a + b)/2M For the Thue-Morse QW structure the 
function f{q,N) in Eq. ([3]) tends to zero when N ^ 00 
as a power of N at any q except certain singular values. 
The latter form a series 



2q = Gh = ^, h^Q,±l. 
d 

with the structure factor given by^i 



(TM), 



COS 



TTah\ 2 f 7r6/i\ 

~2l ) V "2J j 



(10) 



(11) 



The resonant Bragg conditional for both the Fibonacci 
and Thue-Morse structures is formulated as 



c ~ 2 



(12) 



where G stands for Ghh' in the Fibonacci case and for 
Gh in the Thue-Morse case, see Eqs. ([5]) and ((TU]), re- 
spectively. Of course, one can impose a similar con- 
dition for non-singular wavevectors contributing to the 
structure factor of the Thue-Morse sequence. Since in 
this case the value f{q, N) decreases with increasing N 
the corresponding system is far from being an efficient 
exciton-polaritonic structure. This is the reason why 
we do not consider here, e.g., the period-doubling se- 
quence ABAAABAB . . determined by the rule A 
AB, B — > AA, which has no Bragg peaks except for the 
trivial one at q = 0. 
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For the sake of completeness, we also analyze a slightly 
disordered structure with the long-range order main- 
tained and the QW positions defined by 



(13) 



where the deviation Sz^ is randomly distributed and 
defined by the vanishing average, {6zm) = 0, and the 
dispersion = {6z'^). The structure factor f{q) — 
limTv^oo fil, of such a lattice averaged over the dis- 
order realizations has the form 



h 



Gh^2TTh/d. (14) 



The dispersion of f{q,N) tends to zero with N — > oo, 
and Eq. ([M]) provides a good estimation of the struc- 
ture factor for any fixed disorder realization whenever 
N > 10. The long-ranged correlations of QW positions 
are preserved by (|13p. and the Bragg diffraction is pos- 
sible with the same diffraction vectors as in the periodic 
lattice. However, the structure-factor coefficients drop 
drastically with the growth of (T^. The exponential fac- 
tor in is equivalent to the Debye- Waller factor caused 
by the thermal motion of atoms in a crystalline latticci^ 
Since the geometry of MQW structures under study is 
now described and the resonant Bragg condition is im- 
posed we proceed to the optical reflection spectra. 
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III. TWO REGIMES IN OPTICAL 
REFLECTION FROM RESONANT BRAGG 
STRUCTURES 

Two different regimes have been revealed in opti- 
cal reflection from the periodic resonant Bragg QW 
st ructure s. ^"^'^^ For small enough numbers of QWs, N ^ 
a/wo/Fq (superradiant regime), the optical reflectivity 
is described by a Lorentzian with the maximum value 
[iVFo/(iVFo + r)]2 and the halfwidth NTq + F. For a 
large number of wells, N > ^wo/Fq (photonic crystal 
regime), the reflection coefficient is close to unity within 
the exciton-polariton forbidden gap and exhibits an oscil- 
latory behaviour outside the gap. The calculations pre- 
sented below demonstrate the existence of similar two 
regimes for the deterministic aperiodic structures. 

A. Superradiant regime 

The numerical calculation of reflection spectra is car- 
ried out using the standard transfer matrix technique'^^ 
given with more details in Sec. |Vl Figure [T] presents 
the reflectivity Rn{(^) calculated for the light normally 
incident from the left half-space z < upon four dif- 
ferent 50-well structures. All the four, namely, the Fi- 
bonacci, Thue-Morse, periodic and distorted periodic 
structures, are tuned to satisfy the Bragg resonant con- 
dition p2p which can be rewritten as A(a;o) = 2d, where 



FIG. 1: (Color online) Reflection spectra calculated for four 
QW structures, each containing A'^ = 50 wells and tuned to 
the resonant Bragg condition 2d — A(tJo): periodic structure, 
with a = b = d (dashed); Fibonacci chain, with a/b = t 
(solid curve); Thue-Morse sequence with a/b — 3/2 (dotted); 
and weakly disordered periodic with — Ao/20 (dashed- 
and-dotted). Note the break on abscissa axis on the panel 
(a) around uu — uuo- Panels (b), (c) and (d) show the same 
spectra in larger scale of the variable {uj — luq) /T'o- Calculated 
for hFo = 50 fieV, hioo = 1.533 eV, F = O.IFq. 



X{lu) — 2t:c/ {ijjni)). This means that, for the Fibonacci 
structure, the value G in Eq. is set to Ghh' with 
/i = l,/i' = and, therefore, for the four structures 
G = 27r/d and their optical properties can be conve- 
niently compared. 

One can see from Fig. [T] that the condition leads 
to high reflectivity of not only the periodic and qua- 
sicrystalline Fibonacci chains2i but also the Thue-Morse 
and slightly-disordered periodic structures. In the region 
|w — wqI > 20Fo, far enough from the exciton resonance 
frequency, the four spectra have similar Lorentzian wings 
with the halfwidth of the order of iVFo indicating the ex- 
istence of a superradiant exciton-polariton mode. The 
magnitude of the wings is governed by a modulus of the 
structure- factor coefficient, For the chosen struc- 

tures this value runs from — 1 (periodic structure) 
and I/gI = 0.95 (distorted periodic) to |/g| = 0.70 (Fi- 
bonacci) and I/gI = 0.65 (Thue-Morse). The spectral 
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wings in Fig. [T] decline monotonously with decreasing 
In addition it should be mentioned that, for the Fi- 
bonacci QW structure tuned to Guh' with h = l,h' = 1 
and analyzed in Ref. [2l], the structure-factor coefficient is 
|/g| ~ 0-9 and the spectral wings in reflectivity are raised 
as compared with those for the Fibonacci structure tuned 
to Gifi. 

In the frequency region around loq the reflection spec- 
tra from the non-periodic structures show wide dips 
where the reflection coefficient oscillates with the period 
of oscillations decreasing as uj approaches coq. As shown 
below, see also Ref^i, the spectral dip naturally appears 
for a multi-layered deterministic system tuned to a Bragg 
diffraction vector with the structure-factor coefficient fa 
smaller than unity, and it widens as the value of 1 — |/g| 
increases. 

For small values of the exciton nonradiative damping 
rate F (lying beyond experimentally available values), 
the fine structure of optical spectra in the narrow res- 
onance region ranged over few Fq is very intricate, see 
Figs. [T{b)-[Hd). All the considered aperiodic structures 
possess a narrow middle stop-band embracing the exci- 
ton resonance ujq- In particular, for the Fibonacci QW 
structure this stop-band is located between ujq — 0.4Fo 
and uo + O.OFq. The spectral properties in the frequency 
range — wqI Fq for F ^ Fq are discussed in Sec. |V] 
in more details. 

In realistic semiconductor QWs the nonradiative decay 
rate is larger than (or comparable with) Fq, and the ma- 
jority of spectral fine-structure features are smoothed. 
The influence of the nonradiative damping on the reflec- 
tivity from the Fibonacci MQWs is illustrated in Fig. ^ 
Thin curve on the upper panel is the same as that on 
Fig. [1] and calculated for F = O.IFq while the thick curve 
corresponds to more realistic value F = 2Fo. One can 
trace the smoothing of sharp spectral features with in- 
creasing F. However, some of these features may still be 
resolved by means of the differential spectroscopy widely 
used in studies of bulk crystals and low-dimensional 
structures. Panels (b) and (c) of Fig. [3] present the first 
and second derivatives R'{uj) = dR{uo)/dLd and R"{lo), 
respectively. The differential spectra allow one to en- 
hance the spectral peculiarities poorly resolved in the 
ordinary spectrum of Fig. [2{a). 



OS 



1.0 

0.8 
0.6 
0.4 
0.2 
0.0 



■ (a) 


1 1 1 




r-o.ir„ - 


N=50 , I 


r-2r„ _ 















•e 

03 



3 

43 



- — ^ 





















-50 



-25 



25 



50 



(«-co„)/r„ 



FIG. 2; (Color online) Differential reflection spectra cal- 
culated for Fibonacci QW structures containing TV = 50 
wells. Panel (a) shows reflection spectrum R{uj) calculated 
for r = O.lFo (thin curve) and F = 2ro (thick curve). Panels 
(b) and (c) demonstrate the first- and second-order differen- 
tial spectra R' (uo) and R"((-o) in arbitrary units for F = 2ro. 
Other parameters are the same as in Fig. [1] 
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B. Photonic-quasicrystal regime 

The superradiant regime holds up to ~ 100 and is 
followed by a saturation with the further increase in the 
number of QWsi^ In this subsect ion we study very long 
QW chains with N ^ y^wo/Fq. The calculations for 
A^ = 1000 are presented in Fig. [3] for Fibonacci MQWs 
and Fig.[3]for Thue-Morse MQWs. Figure[3]allows a clear 
interpretation of the spectral properties of excitonic po- 
laritons. Two symmetric stop-bands^ ^ standing out be- 
tween numerous sharp maxima and minima are clearly 
seen in the spectra. Figure El^b) shows the spectrum cal- 



FIG. 3: (Color online) Reflection spectra calculated for a Fi- 
bonacci QW structure containing 1000 wells. The parameters 
used are the same as those in Fig. [1] except for the nonradia- 
tive decay which now is F = 0.2ro. Upper and lower panels 
correspond, respectively, to the exact calculation and calcula- 
tion in the two- wave approximation. Vertical arrow at a; = wo 
in panel (a) indicates the narrow reflectivity stop-band, which 
cannot be described by the two-wave approximation. 



culated in the two-wave approximation taking into ac- 
count only three terms in the sum ^ , namely, the terms 
with 2q = ±Gi_o and 2q = 0, see Sec. |IV] for details. 
Comparing Figs. [3Ka) and [S^b) we conclude that a lot 
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of spectral features are reproduced as the interference 
fringes in the approximated spectrum. However, this ap- 
proximation lacks an adequate description of the reflec- 
tion spectrum around the exciton resonance frequency. 
The middle stop-band at w = cjq found in Fig. [ijb) re- 
veals itself also in Fig. ^a) where it is indicated by a 
vertical arrow. 

Thin solid curve in Fig. 2] illustrates the exact re- 
flectivity calculation for the Thue-Morse structure with 
N = 1000. Thick solid curve is calculated in the two- 
wave approximation in the limit iV cxd so that all the 
interference fringes are smoothed due to the finite value 
of r — > -1-0. We have checked for the Thue-Morse se- 
quence that, for N — 1000, the two-wave approximation 
works well outside the narrow interval — wqI Fq. The 
fine features around loq are again beyond the scope of the 
two- wave approximation. The spectrum for 1000 periodic 
QWs, dotted curve in Fig. 31 is presented for comparison 
in order to emphasize similarities and differences between 
the optical spectra of periodic and aperiodic systems un- 
der consideration. 

Figures [TJ3] form a computational data base for the 
physical interpretation of the reflection spectral shapes. 
This can be done in terms of the two- wave approximation 
(Sec. IIVP and scaling and self-similarity (Sec. |V|. 
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FIG. 4: (Color online) Reflection spectra calculated for Thue- 
Morse and periodic QW structures. Thin solid curve corre- 
sponds to the exact calculation for = 1000 Thue-Morse 
QWs, thick solid curve is calculated in the two-wave approx- 
imation for A'^ oo. Dotted curve represents the A'^ = 1000 
periodic MQWs with the period d = A(aJo)/2. The parame- 
ters used are the same as those in Fig. [S] 



IV. OPTICAL SPECTRA IN THE TWO- WAVE 
APPROXIMATION 

In this section we apply the two-wave approximation 
for derivation of the exciton-polariton dispersion and the 
reflectivity spectra of the aperiodic resonant Bragg MQW 



structures. The derivation is performed for the particular 
case of a canonic Fibonacci chain but the results can be 
straightforwardly generalized to noncanonical quasicrys- 
tal sequences and other deterministic aperiodic multi- 
layered structures. 

The electric field of the light wave propagating in the 
MQW structure satisfies the following wave equation^! 



E{z) 



2gFo 



iF 



^(5(z - Z,n)E{Zyn) ■ 
m 

(15) 

Assuming the resonant Fibonacci MQWs to contain a 
sufficiently large number N of wells we replace the struc- 
ture factor f{q,N) by its limit f{q) given by Eq. H]). 
This allows one to approximate solutions of Eq. (jlSp as 
a superposition of the "Bloch-like" waves 

Ek{z) = e'(^^'=--')^i?K-G,„ . (16) 

h,h' 

As distinct from the true Bloch functions in a periodic 
system, here the countable set of vectors Ghh' is enumer- 
ated not by one but by two integer numbers, h and /i', 
see Eq. ([5]). Substituting Ek{z) into Eq. (fTS]) . multiply- 
ing each term by eyi'p[—i{K — Ghh')z\ with particular h, h' 
and integrating over z we obtain 



hh' 



'2qi]EK-G,,, 



99' 



where 



e(^) 



To 



(17) 



(18) 



diujQ — Lo ~ iF) 

Note that, throughout this paper, we focus on a fre- 
quency region |a; — ljoI *C around the exciton reso- 
nance. 

In accordance with Eq. (fT^ we consider the Fibonacci 
QW structure tuned to the Bragg resonance 

^0 _ Ghh' 



-Tib = qo 



(19) 



In the two-wave approximation, only two space harmon- 
ics K and K' — K ~ Ghh' = K — 2qQ are taken into 
account in the superposition p^ . A necessary but not 
sufficient condition for validity of this approximation is 
the inequality 



\qo -Kl-^qo 



(20) 



Then the infinite set (fT7|) is reduced to a system of two 
coupled equations 

{q-K + OEK+CfL'EK'=0, (21) 
^hh'EK + {q + K- 2qo + OEk' - . 

The two eigenvalues K corresponding to the frequency uj 
are given by 

= go ± Q , Q = + g - q^Y - e\fhh'V ■ (22) 
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The criterion (12011 is then rewritten in the form 



max{|ti;o - T} > Tq 



h+h'/T 



(23) 



At fhh' ^ 1 this dispersion equation reduces to that 
for the periodic resonant Bragg MQWsi^ 



Q = qo] 



where 



UJ — LUq 
. ^0 



A 



o; — 0^0 + ir 



2roWo 



The edges, and cjj^, of two symmetrical band 

gaps in the Fibonacci QW structure are obtained from 
Eq. ^ by setting K^^'^ — qo or, equivalently, Q — 0. 
The result reads^i 



± _l_ A / + Ifhll' 



2{h + h'/r) ' 



(24) 



UJ. = cjQ ± A. 



1-1// 



hh' 



2{h + h'/T) 



As shown in Appendix, for the exciton-polariton waves 
at the bandgap edges located at the point K ~ qo, a.n ad- 
mixture of other space harmonics has no remarkable in- 
fluence and these edges are well-defined for the resonant 
QW Fibonacci chains. For each frequency lying inside 
the interval between the edges cJj" and Wjjj, the two- wave 
approximation gives two linearly independent solutions 
^ qo±Q with nonzero Q. An exciton-polariton 
wave induced by the initial incoming light wave is a su- 
perposition of these two Bloch-like solutions. They can 
be coupled by the diffraction wavevector Gggi satisfying 
the condition 2Q = Ggg' . If the corresponding structure- 
factor coefficient fgg> is remarkable one should include 
into consideration mixing of the waves which com- 
plicates the comparatively simple description of exciton 
polaritons. In the approximate approach we will ignore 
the diffraction-induced mixing between the waves 
and analyze the validity of this description by comparing 
the exact and two- wave calculations. 

In order to derive an analytical expression for the re- 
flection coefficient from an A^-well chain sandwiched be- 
tween the semiinfinite barriers (material B) we write the 
field in the three regions, the left barrier, the MQWs and 
the right barrier, as follows 



^_e-"5^(e'9o2 -I- <-_e^'«o^) 



(z < 0) , 

(0 < z < Nd) 
{Nd < z) . 

(25) 



Here rjv and tj\r are the reflection and transmission am- 
plitude coefflcients, E± are the amplitudes of the "Bloch- 
like" solutions, and 



c± 



hh' 



q-qo±Q + ^ 



Values of rN,tN, E- are related by imposing the 
boundary conditions which are continuity of the elec- 
tric field E{z) and its first derivative dE{z)/dz at the 
points zi = and z^r. If the number of wells N co- 
incides with N = Fj + 1, where Fj is one of the Fi- 
bonacci numbers, then the product G^h'^N differs from 
an integer multiply of 27r by a negligibly small value, 
Sn = —Ghh's{l — tY^^t^'^. In this case the phase factor 
exp (iGhh'ZN) can be replaced by unity and the straight- 
forward derivation results in surprisingly simple expres- 
sions for the reflection coefficient, 



TN = 



hh' 



qo - q- C-iQcot (QNd) 



(26) 



and for the transmission coefflcient, 
tN - -irArO/K/w,' siniQNd)]. 

Numerical calculation demonstrates that in the region 
(|23p the two-wave approximation is valid even for the 
Fibonacci structures with 7^ Fj + 1 provided that 

> 20 and the mesoscopic effects are reduced. More- 
over, Eq. (|26p for the reflection coefflcient can be applied 
to other deterministic aperiodic systems including the 
Thue-Morse and weakly disordered periodic structures. 
It suffices for the structure to be characterized by a sin- 
gle value of the structure factor at particular vector which 
can play the role of the Bragg diffraction vector. In the 
crystalline case where fhh' = 1 analogous approximation 
describes the nuclear resonant scattering of 7-raysi^ 

The two-wave approximation can be generalized to al- 
low for the dielectric contrast, i.e., the difference between 
the dielectric constant of the barrier, n^, and the back- 
ground dielectric constant of the QW, n^. In a structure 
with Ha 7^ rih the stop-band exists even neglecting the ex- 
citon effect, Fq = 0. The excitonic resonance leads to the 
splitting of this single stop-band to two ones. In the peri- 
odic case the highest reflectivity is reached when the two 
stop-bands touch each other and form a single exciton- 
polariton gap. This is the effective Bragg condition for 
the periodic structure with the dielectric contrast. In 
the Fibonacci case when \fhh'\ < 1 the stop bands never 
touch each other and the Bragg condition means that the 
sum of their widths reaches a maximum. For the real- 
istic case of a small contrast, \na — nh\ <Si na,nb, this 
condition is equivalent to the tuning of the exciton reso- 
nance frequency luo to one of the edges of the stop-band 
found at Fq = 0, similarly to the condition for the peri- 
odic structures] ^^1**° . Note that the reflectivity spectrum 
taken from the Bragg MQW structure with the dielectric 
contrast is always asymmetric. 

For periodic resonant Bragg MQWs, in the superradi- 
ant regime N <^ ^^luo/To, Eq. is readily transformed 
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to the well-known resulljiS 



iNTo 



ojQ-Lj- i{NTa + r) 



(27) 



The pole sA uj = ujq — i(iVro + F) is the eigenfrequency of 
the superradiant mode. In general, the eigenfrequencies 
w^'^ of a MQW structure are represented by zeros of the 
denominator in Eq. (|26p . Since the structure is open 
the eigenfrequencies are complex even in the absence of 
nonradiative damping, F = 0. Values of lo^^^ lying in the 
region 



\U} - UJq\ < 



hh' 



y/h + h'/r 



(28) 



but outside the narrow interval (|23p can be easily found 
taking into account that, in this region, the difference 
qo — q in Eqs. ((22)) . ([26|l can be neglected as compared 



with ^ so that one has Q = ^-y/l — Ifhh'l"^ and 



ui 



(0 



loq — iT — liVFo 



arctanh(^l - lAvp) + iirl 
/ = 0,±1,±2... (29) 



Equation 



determines at I = the frequency of the 

fhh' 

this frequency is given by 



superradiant mode. In the particular case 1— <C 1, 



(0) 



: ujo ~ iNTo 



hh' 



(30) 



-2 - 



-3 -. 
o 

-40'- 
-45 ^ 



-10 



10 



Re ((a-(o^)/r^ 



FIG. 5: (Color online) Complex eigenfrequencies of exciton 
polaritons in Fibonacci 56-QW structure calculated exactly 
(filled symbols) and in the two-wave approximation (empty 
symbols) . Note the break at the ordinate axis. Vertical arrows 
indicate the edges of inner stop-band. Calculated for F = 
and other parameters same as in Fig. \T\ 

Figure \5\ shows the eigenfrequencies of a 56-QW Fi- 
bonacci structure calculated exactly (filled symbols) and 



from Eq. (|29|) (closed symbols). The exact calculation 
is performed using the system of coupled equation for 
excitonic polarization P„i in the quantum well^ 



N 



{ojQ-uj- iF)P„i -I- ^ 



>IP„,,, m=l...N. 



(31) 

One can see from Fig.[5]that the two- wave approximation 
excellently describes the superradiant mode as well as 
some of the subradiant modes lying far from ujq on the 
complex plane. The approximation breaks in the region 
close to ajQ where more sophisticated analysis is required, 
as presented in the next Section. 



V. OPTICAL SPECTRA NEAR THE EXCITON 
RESONANCE FREQUENCY 

A. Trace map technique for Fibonacci and 
Thue-Morse quantum well structures 

In this section we concentrate on the narrow frequency 
region |w — wqI ~ Fq. A convenient and powerful tool for 
such study is the transfer matrix method. The substitu- 
tion rules ([1]) lead to recursive relations for the transfer 
matrices providing all the essential information about the 
spectral properties of MQWs. 

In the following we use the notation Tj for the Fi- 
bonacci chain containing = Fj QWs starting from the 
trivial chain J-'i that consists of one segment A. Sim- 
ilarly, TAij is the Thue-Morse sequence with N — 2^ 
QWs starting from TMq = A. The transfer matrix 
Mj through the whole structure or TAij is given by 
a product of the matrices Mqw, Ma, Mb of transfer 
through a QW and a barrier of length a or b, respec- 
tively, with the order established by the chain definition 
([1]). In the basis of electric field E{z) and its derivative 
—q^^dE{z)/dz the transfer matrices are as follows^ 



QW 



V 



f cos qd 
\SVL\qd 



(1 
\2S 1 



— sin qd 
cos qd 



Fo 



LjQ — uj — iT 



X> = A,B;d 



(32) 



(33) 



We will here restrict ourselves to the limit of zero nonra- 
diative decay, F = 0, in which case the transfer matrices 
are real. The transmission and reflection spectra, Tj{uj) 
and Rj{uj), are given by^ 



T,(c.)-l-i?,(^) 



1 



(34) 



Here the quantities Xj and yj stand for the half-trace 
{Mj^ii+Mj,22)/2 and half-antitrace {Mj,2i- Mj^i2)/2 of 
the matrix Mj, respectively. 
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In order to reveal the behavior of exciton polaritons 
in aperiodic MQWs it is instructive to calculate the po- 
lariton dispersion in the approximants^^ of the aperiodic 
chains containing the periodically repeating sequences Tj 
or TM.j. In such periodic systems the polariton band 
structure consists of allowed minibands and forbidden 
gaps. The gaps are found from the conditional 

|x,M|>l. (35) 

To proceed to the analysis of the pattern of allowed and 
forbidden bands we note that the half-traces Xj of the 
substitution sequences satisfy closed recurrence relations, 
also termed as trace maps.— For the Fibonacci and Thue- 
Morse chains, the trace maps read 

Xj+i = 2xjXj-i — Xj^2 (Fibonacci) , (36a) 
Xj+i = 4a;^2_i(xj- - 1) + 1 (Thue-Morse) . (36b) 

Consequently, the polariton energy spectrum is deter- 
mined by the general properties of nonlinear transfor- 
mations (I36p and the initial conditions specific for the 
QW transfer matrices ((32|) . The trace maps are effec- 
tive for studies of the spectral properties of deterministic 
aperiodic structuresi^i^i 
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(a) 



0.92 0.94 0.96 0.98 1 1.02 1.04 



FIG. 6: (Color online) (a). Exciton-polariton allowed (thin 
lines) and forbidden (thick stripes) bands in periodically re- 
peated Fibonacci sequences of the order j — 1 ... 13. (b) 
Bands for j = 11 ad j = 13 in the spectral range around 
the frequency uj — loq. (c),(d) Bands for j = 11 and 
j = 13, respectively, in a large scales near the frequency 
UJ = ujo + 0.935ro indicated by the vertical line. Calculated 
for r = and other parameters same as in Fig. [T] 



B. Scaling of band structure and transmission 
spectra in Fibonacci structures 

We will now analyze the band structure and the 
transmission spectra in Fibonacci QW structures. For 
the Fibonacci lattices the trace map (|36ap possesses an 
invariant^ 



values at particular frequencies lUs where the reflectivity 



I — Xj + Xjj^i + Xjj^2 



2XjXj+iXj+2 — 1 ■ 

In the QW structure this invariant can be reduced to 

I{uj)^ S^{uj)sm^[qb{T-l)]. (37) 

The resonant behavior of this invariant as a function 
of frequency indicates that the band structure for Fi- 
bonacci QW chains may be very complex in the region 
|w — wol ~ Fq. The band calculations are presented in 
Fig. IHa), where the black stripes and horizontal lines 
show, respectively, the forbidden and allowed bands for 
different values of the structure order j ranging from 
j = 1 to j = 13. Figures[nKb)-[nKd) represent this band se- 
quence for j = 11 and j — 13 in different frequency scales. 
Panel (a) demonstrates that two broad band gaps are al- 
ready present for 21 QWs (j = 8). With increasing j 
their edges very quickly converge to the analytical values 
((24)) shown by the gray rectangles in Fig. [Bfa). A narrow 
permanent middle band gap at — 0.4ro ^ lu — luq ^ O.OFq 
is well resolved in the scale of Fig. (Hb). 

The other forbidden bands depicted in Fig. [S] can be 
interpreted in terms of two formation mechanisms. The 
first mechanism is related to the two-wave approxima- 
tion. In this approximation the half-trace of the trans- 
fer matrix Xj reaches minimum (—1) or maximum (-1-1) 



,(2w) 
N 



vanishes. One can check from Eq. (pS)) that at these 
frequencies the product QNd is an integer number of tt. 
In the vicinity of one can approximate x*-^"-* {uj) by 
ns[l — Us{uj — Ws)^], where = ±1, and Us is a positive 
coefficient. Near ojs the exact function Xj (w) = Tr(Mj)/2 
differs from x'^'^'™\uj) by the correction 5x{uj) which can 
be approximated by ns[cs + Vs{uJ — uJs)] where Cs, Vs are 
additional constants. As a result the behavior of the half- 
trace can be presented in the form 



x{uj) 



4ms 



2u, 



If Cs + (Vg/Aus) is positive then the periodic system has 
a gap at uj'g — Us + {vs/2us). 

The second mechanism of gap formation is related 
to localized exciton-polariton states rather than to the 
Fabry-Perot interference and can be treated in terms of 
self-similarity effects. Particularly, in the frequency range 

— (jJqI ~ Fq the number of stop-bands increases while 
their widths tend to zero as ^ oo. As a result, the 
sequence of the allowed and forbidden bands becomes 
quite intricate, see Figs.[SKb)-lSl[d), and locally resembles 
the Cantor set 4^ The most striking result in Fig. ^h) is 
similarity of the band structure of the approximants with 
j = 11 and j — 13. On the other hand, the spectrum for 
j = 13 has a lot of narrow band gaps not resolved in 
the scale of Fig. [6Kb) . Figures [IJc) and ^d) present the 
same spectra in larger scales near the right edge of the 
middle band gap, with the scale for j = 13 being A+ w 8 
times larger than that for j = 13. Matching the bandgap 
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positions we prove an existence of the spectral scaling 
in the Fibonacci QW structures. The scaling index A+ 
specifies the ratio of the widths of spectral features of 
the structures with the order differing by two. The scal- 
ing properties hold not only for the band positions but 
for the whole curves Xj{uj), as Fig. [7] demonstrates for 
j = 11, 13 and 15. One can see that the curves plotted 
in the proper scales repeat each other almost exactly. 

Now we turn from the transfer matrix traces anal- 
ysis, to the complex eigenfrequencies ([31]) and to the 
transmission spectra Tj{uj). The spectra Tj{u}) for the 
resonant Fibonacci structures are shown in Fig. [5] for 
j = 11 (iV = 89 QWs) and j = 13 (N = 234 QWs) 
in the frequency region adjacent to the right edge of the 
middle band gap. The real parts of the complex eigenfre- 
quencies are indicated by the vertical lines. The abscissa 
scales on Figs.[8l^a) and[8l^c) are the same as in Figs.|6l^d) 
and[6l^c). Examining Fig.[8]we conclude, that the scaling 
properties revealed by Figs. [6] and El are also manifested 
in the optical spectra. Indeed, comparing Fig. [5]^ a) and 
Fig. El^b) we find that the general shape of the spectra 
remains similar although more details appear when N 
grows. On the other hand, the relative distances between 
the transmission peaks for j — 13 in the large scale agree 
with those for = 11 in a smaller scale, cf. Fig. [H^a) 
and Fig. [Sjc) . The positions of the real parts of complex 
eigenfrequencies correspond to the peaks in the trans- 
mission spectra and exhibit the same scaling behaviour. 
Such behavior is also demonstrated at the left edge of the 
middle bandgap w — wq ~ — O.STFq, it is characterized by 
the scaling index A_ ^ 16. 

The self-similarity of band structure of Fibonacci se- 
quences with the orders j differing by 2 can be related 
to the so-called "band-edge" cycle of the trace map 
P6ap 4^ This can be done by the following considera- 
tion. If the half-traces Xj for three successive orders 
j = jo, jo + 1, jo + 2 are interrelated by 



quencies = LUj:'{jo oo) = Lu'£^{jo 
values ujI « 0.94ro and uj*_ w -0.37ro . 



2x 



, Xjg+2 



(38) 



Jo 



then according to the recurrent equations (|36ap the val- 
ues of Xj for j > jo form a periodic sequence 



oo), with the 



I 




0.935 0.936 0.937 0.938 0.939 0.940 0.941 
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0.9350 0.9352 0.9354 0.9356 

{uj - ujQ)/ro 



2x* 



2x* 



2x* 



(39) 



where x* = Xjg. Although the sequence (139]) repeats 
itself only after 6 iterations, the length of the cycle for 
the absolute value of the trace is two. 

Considering Xj as functions of the frequency uj we in- 
troduce solutions LlJ^°■^jo) and w'^^'(jo) of the first and 
second equations (|38p . The numerical calculation shows 
that, for each Jq, in the vicinity of luq there exist two pairs 
of solutions, uj'^\jo), and Lu'^\jo)- Moreover, values of 

uj'f\jo) and w+^(jo) or wL°^(jo) and wL^^(jo) merge with 
increasing jo, and one can introduce the asymptotic fre- 



FIG. 7: (Color online) Trace scaling for the Fibonacci struc- 
tures. The panels (a), (b) and (c) show the half-trace Xj{uj) 
for j — 11, 13, 15, respectively. The filled ribbons indicate the 
regions of polariton band gaps, where |j:j(cj)| > 1. Calculated 
for r = and other parameters same as in Fig. [T] 



Our analysis made for Fibonacci structures with dif- 
ferent pair values of {h, h') shows that they also demon- 
strate analogous scaling behaviour in the vicinity of wg- 
The distance \uj\_ — uli*_ \ between the scaling frequencies 
decreases with the growth of the barrier thicknesses. The 
scaling indices increase, when the middle band gap be- 
comes narrower. We have established that there exists 
the following equation relating the scaling indices and 
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(ffl-«„)/r„ 

FIG. 8: (Color online) Transmission spectra of the Fibonacci 
quantum-well sequences of the order j =ll(a), and j =13 
(b,c), containing 89 and 233 QWs, respectively. Vertical 
lines indicate the real parts of complex eigenfrequencies of 
the structures. Calculated for F = and other parameters 
same as in Fig. [T] 



ing features in the region |tj — wqI ~ Fq. At very high 
values of > 1000 the optical spectra become intricate 
even at jtj — tJol 3> Fq. However, in the stopband re- 
gions (tj~,j,a;r) and {ui^, uJ'^^^^) the reflection coefficient 
Rn{uj) — kwCw)!^ remains close to unity. The spectral 
pattern within the interval (wj^, Wj^J strongly depends on 
values of T and N. For any nonradiative damping rate F 
there exists a finite value of the number of wells, N{T), 
which separates the structures into two categories. For 
N > A^(F), the reflection spectrum is independent of N, 
i?Ar(a;;F) « i?oo(w;F), and determined by the exciton- 
polaritons localized within the area < z < 2;Ar(r) and 
sensitive to the initial light. For QW structures with 
N < N{T), the light wave reaches the back edge of the 
structure, reflects from this edge, propagates back and 
participates in the Fabry-Perot interference resulting in 
the oscillating reflectivity. This regime, N < A(F), is 
well described by the two- wave approximation except for 
the narrow region near the exciton resonance frequency 
where the condition ([^5]) is not satisfied. With decreas- 
ing F the critical number A'^(F) infinitely increases while 
the spectrum i?oo(aj;F) continuously varies as F ^ -1-0 
and shows no saturation behaviour. 



the frequencies ujj.: 

A± « A1+A2 J/(c^|) = Ai+ , , \sm[q{u;o)b{T-l)] 

V \Uo—UJ^\ 

(40) 

The coefficients Ai and A2 are found to be close to 3 
and 4, respectively, and are independent of the structure 
parameters. Since the value of \uj*l — ljo\/To w 0.94 ex- 
ceeds |w* — wo|/ro « 0.37 the scaling coefficient A+ « 8 
is smaller than A_ « 16. 

We have also analyzed the spatial structure of 
the excitonic polarization of the eigenstates sat- 
isfying Eq. (|3ip . In particular, this distribution 
is characterized by the participation ratio p — 

I]m=i l-P™l'*/(Em=i l-PmP) : where the sum runs over 
QW-lattice sites4^ The parameter p is a measure of the 
state localization-delocalization: for a completely delo- 
calized state p « 1/A^, whereas p ~ 1 for a state tied 
to a single site. In the periodic Bragg structure with 
d = X{luq)/2, the superradiant mode is described by the 
eigenvector P,„ — {—1)™'Pq and the participation ratio 
p — 1/N. In the Fibonacci QW chain the participation 
ratio remains small for the superradiant mode as well 
as for the subradiant modes described by the two-wave 
approximation (|29p . On the other hand, there exist local- 
ized states with large values of p as well as intermediate 
states. The eigenstates with strongly localized charac- 
ter belong to the region \uj — loqI covering few Fq. The 
excitonic polarization of such states is concentrated 
on a small fraction of the QW chain and has a complex 
self-similar structure.^'* 

Up to now we have limited ourselves only to the scal- 



C. Transmission spectra of the Thue-Morse 
quantum well structures 

We now turn to the Thue-Morse structures. The po- 
lariton band-structure calculations performed for this 
system lead to qualitatively similar conclusions: two- 
wave band gaps are already formed for small j , a middle 
narrow band gap is always present, and a complicated se- 
quence of the allowed and forbidden bands arises around 
luq. However, the Thue-Morse structures have a very in- 
teresting specific properties, most brightly manifested in 
their transmission spectra. 

The transmission spectra are presented in Fig. [5] for 
different orders j changing from 4 (A^ ~ 16 QWs) to 
7 (A^ = 128). The spectrum has a complex structure 
with narrow peaks even for N = 16, see Fig. HJa). Fig- 
ures|ni[b)-[Hd) show evolution of the spectra with increas- 
ing the number of QWs. The trace map (|36bp alone is 
not sufficient to obtain transmission spectra. Therefore, 
to analyze the spectra we use the standard properties of 
the trace map (j36b|) and the antitrace maps^ 

yj+i = 2xj^i[{2xj - l)yj-i + yj-i] , (41) 
Vj+i = 2a;j-i[(2a;j - l)yj-i + yj-i] . 

Here the half-antitrace j/j corresponds to the structure 

TMj obtained from TAij by the barrier interchange 

A^B, e.g., TM2 = ABBA and TAla = BAAB. Con- 
trary to the Fibonacci case, the trace map (j36b[) for the 
Thue-Morse structures has no cycles of type (P^ . In- 
stead of such cycles, Eqs. (j36bp and (HT1) have the follow- 
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FIG. 9: (Color online) Fine structure of the transmission spec- 
tra calculated for Thue-Morse quantum-well sequences of or- 
der j = 4, 5, 6, 7 (AT = 16 . . . 128). Calculated for F = and 
other parameters the same as in Fig. [T] 



ing propertji^ 



if Xj^2{^) = then 



Xj{uj) = a;j_|_i(ti>) = . . . = 1 
Vji'^) = Vj+ii^) = • . . = . 



(42) 

As a consequence, the structure becomes transparent at 
this particular frequency: Tj{uj) = Tj^i{uj) = ... = 1. 
The positions of the spectral features follow then from 
Eq. (|^^ . For example, the edges of the inner band gap for 
the Thue-Morse structure, ujl = ojq — CSSFo and loh = 
Ci;o + 0.47ro, can be found from the equations xi{bJL) = 
and X2{ojr) = 0, respectively. 

Whenever the half-trace Xj-2{'.^) vanishes at some fre- 
quency w = wi and therefore, Xj{uji) = 1, there exist two 
neighboring frequencies, uj2 < oJi and uj^ > uji, such as 
Xj{t02.,3) = 0. As a consequence, Xj+2{^2,3) = 1 and the 
transmission coefficient for the structure TA4j+2 at these 
frequencies reaches unity. Thus the number of trans- 
mission peaks increases with the growth of the structure 
order: each unitary peak in the spectrum of the struc- 
ture TAij (i) persists in the spectra of the structures 
TA4j+i,TA4j+2, ■ ■ ■ of higher orders and (ii) leads to 
appearance of two more adjacent peaks for the structure 
TA4j+2- An example of such "tree" of trifucations is 
indicated in Figs. [9l^b)-[9l[e) by dashed lines. The sin- 
gle transmission peak of the structure TM4 {N = 16) 
a,t uj ojQ + 2.IF0 has generated the two peaks at 
uj fa wo + l.SFo and uj « cjq+Z.GFo for the structure TMe 
{N = 64), cf. Fig.HJb) and Fig.[ltd). The characteristic 
widths of the spectral features tend to zero proportional 
to a power the structure length N"" = 2"-'", where a 
is a frequency-dependent positive index. Transmission 



spectra can also have non-unity maxima, see Fig. ^e). 
These peaks do not correspond to any special values of 
Xj and Uj and their positions depend on j. The spatial 
distribution of the electric field on the frequencies with 
unitary transmission has a so-called "lattice-like" shape, 
specific for the Thue-Morse structures4^ 

The spectra presented in Fig. [9] are calculated for the 
Bragg structure with q{u!o){a + b) — 27r. An interesting 
property of the Thue-Morse lattices, satisfying the Bragg 
condition q(LUQ){a+b) = tt, 27r . . ., is the mirror symmetry 
between the transmission spectra T(w) and T{uj), of the 
structures TAij and TAij, holding in the region \uj — 
LJ{)\ <JC ^/Tqujq in which case q 



q{uj) in Eq. can 
be approximated by (7(^0)- The spectra are symmetrical 
with respect to uj ^ loq: T[ijj) ~ T(2luo — uj). 



VI. CONCLUSIONS 

We have developed a theory of exciton-polaritons in 
photonic-quasicrystalline and aperiodic MQW structures 
having two different values of the interwell distances, A 
and B. The approach based on the two-wave approxi- 
mation has been extended to derive both the dispersion 
equation and analytical formulas for the reflectance and 
transmittance spectra from aperiodic MQW sequences. 
This approximation successfully describes the pattern of 
the optical spectra including the pair of stop-bands and 
interference fringes between them but stops being valid in 
a spectral range |u; — wqI Eq near the exciton resonant 
frequency ujq. In the Fibonacci QW chains with small val- 
ues of the exciton nonradiative damping rate, E < Eq, the 
transmittance spectra and polariton band structure re- 
veal, at two edges of the narrow inner band gap, a compli- 
cated and rich structure demonstrating scaling invariance 
and self-similarity. It has been shown that this structure 
can be related to the "band-edge" cycle of the trace map. 
The fine structure of optical spectra in the Thue-Morse 
MQWs can be interpreted in terms of zero reflection fre- 
quencies, or frequencies of unity transmittance4i 



Acknowledgments 

This work was supported by RFBR and the "Dynasty" 
Foundation - ICFPM. The authors thank G. Khitrova 
and H.M. Gibbs for helpful discussions. 



Appendix 

Here we will analyze the second order of the perturba- 
tion theory at the point K ~ —K' = G'^^,/2 = go and 
confirm the high accuracy of the stop-band edges defined 
by Eqs. (El. 

We consider the terms in Eqs. pT]) proportional to the 
structure-factor coefficients //t_g h'-g' ^ ^ perturbation. 
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Then keeping the second order contributions we obtain 
the following equations for the amplitudes and Efci = 



+ xd'C22)]E- 

(Al) 



Xihh' + xd^(2i)EK + [q^ -K^+ x(l + xd'<:22)]E-K = 



Here x^M^ 2qTo/[d{LUo - to + iF)], 

A _ y- ^ ^ l/gg' P 

(5,3')5^(0,0) 



C2I - CI2 



fh-g,h'-g' fgg' 

(S,S')7^(0,0) 



(A3) 



We remind that, for the Fibonacci chains, the structure- 
factor coefficients are given by 



_ sin5'gg> 
Jgg' ^ Q 

•^99' 



T — 2 \ ITT 

exp(i— — S-ggJ, Sgg^ = -^^(rg'-g). 

(A4) 

The applied perturbation theory differs from the stan- 
dard one by the presence of terms in the corresponding 
sums with the denominators {K— Ggg')^ — arbitrarily 
close to zero. However the sums are finite because the 
smallness of the denominator at particular values of g and 
g' is compensated by much smaller values of the numer- 
ators for the same values of g,g'. The convergence of the 
sum (|A2p for pairs g, g' where Gggi fa or Gggi ~ Ghh' 
is checked as follows. Taking into account the symmetry 
property \ fgg'\ = \ f-g,-g'\ we can perform the following 
replacement in Eq. (jA2p 
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I/. 



{K-Ggg,f-K^ G^ 



99 



^hh' 



Therefore, this sum converges for the pairs g' with Gggi 
tending to zero. Now let us consider the sequence of pairs 



(g, g') = (h + k,h' + k'), where 

k = mFj, k' — —mFjj^i , 

Fj are the Fibonacci numbers and m is any integer 
±1,±2... different from 0. Taking into account that 



^ - (1 - ry _ ^ (1 - ry-' 

^ 2t - 1 ' r ^ ^ T 



the chosen sequence possesses the properties 

277 (1 - 



G„n' — Ghh' ~ —m- 



SL' {3 -> 00) 



d T 
2 



^20-1) 



2t- 1, 

and, hence, with increasing j one has 

\Ggg. -GhhASlg. « 

which means that the sum 



|Ggg' - Ghh'\ ^Sgg, 



converges. The convergence in Eq. (|A3|) is checked in a 
similar way. 

Numerical calculation performed for the Fibonacci 
QW chain with {h, h') = (1,0) shows that \C,ii\ and \C,2i\ 
are both smaller than 0.1. The prefactor in Eq. (jVip 
is small as compared with fhh' whenever 



, T SttFo /, h' 
max{ LJo - ,r} > — — r ft H 

\fhh'\ V T 



(A5) 



We conclude that, for the wave with K = qo, the 
contributions from the wavevectors Gggi different from 
and Ghh' are negligible within the applicability range of 
Eq. 
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